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Abstract 

It is well known that the zeros of orthogonal polynomials interlace. 
In this paper we study the case of multiple orthogonal polynomials. 
We recall known results and some recursion relations for multiple or- 
thogonal polynomials. Our main result gives a sufficient condition, 
based on the coefficients in the recurrence relations, for the interlac- 
ing of the zeros of neighboring multiple orthogonal polynomials. We 
give several examples illustrating our result. 



1 Preliminaries 

An important and useful relation between zeros of consecutive polynomials 
is the following separation of zeros 

yi < xi < y2 < . ■ ■ < Xn < Vn+l (1) 

where xi, . . . ,Xn and yi, . . . , yn+i are the zeros of the polynomials p„ and 
Pn+i respectively. The most important consequence of this property is that 
it simplifies estimates of the ratio of polynomials, needed for proving the 
asymptotic behavior. Among other things, the interlacing property of zeros 
guarantees convergence of the approximation of zeros by zeros of certain 
special functions using a fixed point method (see |^) and it is equivalent to 
the positivity of weights in quadrature formulas (e.g., [6]). 

The zeros of orthogonal polynomials interlace as a consequence of the 
Christoffel-Darboux formula or the recurrence relation (e.g., [8]). It holds 
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even more general for a Sturm sequence of polynomials. We are interested 
in finding the interlacing property for multiple orthogonal polynomials from 
their recurrence relations. 

In the first subsections 11.11 - 11.31 we briefly introduce multiple orthogonal 
polynomials and their recursion. In Section |2] we summarize known results, 
which are for the case of an AT system, a Nikishin system and a mixed type 
system with two Nikishin systems. We present our main result in Subsection 
12. H which states that the positivity of the coefficients a^i^j in the nearest 
neighbor recurrence relations ensures the interlacing of zeros of neighboring 
multiple orthogonal polynomials. It is worth noting that this condition does 
not depend directly on whether the measures (yU-i, . . . ,/ir-), with respect to 
which the polynomials satisfy orthogonality conditions, form a certain special 
system. Section |3] contains several examples of our result. It gives examples 
of multiple Hermite polynomials (see 13. ip . multiple Charlier polynomials (see 
13. 2p . both kinds of multiple Meixner polynomials (see 13.31 and 13. 7p . multi- 
ple Krawtchouk polynomials (see 13. 4p . and both kinds of multiple Laguerre 
polynomials (see 13.51 and 13. 6p . 

In the following we will use the standard notation, i.e., by n we denote 
the multi-index (ni, . . . ,72^), where rij G N. By \n\ we denote the length of 
the multi-index, i.e., \n\ = J2l=i''^i- 

1.1 Multiple orthogonal polynomials 

As is well known, there are two types of multiple orthogonal polynomials, 
type I and type II (e.g., IH [9]). In this paper we will mainly investigate mul- 
tiple orthogonal polynomials of type II, namely the unique monic polynomial 
of degree \n\ satisfying the orthogonality conditions 

x^Pfi{x) d^j{x) = 0, for < A; < rij - 1, and j = 1, . . . r, (2) 

with respect to r different positive Borel measures (/xi, . . . that are ab- 
solutely continuous with respect to a measure d/i, i.e., for all j = 1, . . . r we 
have dfij{x) = Wj{x) dfi{x). We say that the multi-index n is normal when 
this monic polynomial Pft is unique. 

Type I multiple orthogonal polynomials are polynomials (A^ i, . . . , ,.), 
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where Af^j has degree < rij — 1 such that 

/ x''Afij{x) dnj{x) =0, for < A; < |n| - 2, and j = 1, . . . r, (3) 

and 

^ / dfi,{x) = 1, for J = 1, . . . r. (4) 

These polynomials are uniquely defined if and only if the multi-index n is 
normal. We will use the notation Qft{x) = Y7j=i ^n,j{x)wj{x). 

Multiple orthogonal polynomials are intimately related to Hermite-Pade 
approximants and often they are also called Hermite-Pade polynomials. 

1.2 Recurrence relations for multiple orthogonal poly- 
nomials 

Recall that orthogonal polynomials satisfy a three term recurrence relation 

xpn{x) = a„+ip„+i(x) + hnPn{x) + a„p„_i(x), n > 0, 

with initial values = 1 and p_i = (e.g., [8]). Multiple orthogonal poly- 
nomials also satisfy finite order recurrence relations. There are two types 
of recurrence relations for multiple orthogonal polynomials. Often only the 
multiple orthogonal polynomials are considered with the following sequences 
of indices 

diagonal: n = [n, . . . ,n), 1^1 = ™; 

step-line: n = (m+1, . . . , m+1, m, . . . , m), \n\ = rm + s,0<s<r — 1. 

Any multi-index n on the step-line may be identified by the value of 
the length of n. Observe that every natural number n may be written as 
n = mr + s with < s < r — 1 and then the corresponding multi-index n 
is the one above where the value m + 1 is repeated s times. If all multi- 
indices are normal, then as a consequence of the orthogonality conditions we 
have the following finite order recurrence relation for the polynomials on the 
step-line: 

xPn{x) = Pn+l{x) + an,nPn{x) + a„, 

-|- . . . -|- (ln,n—rPn—ri^X^y ^ 0, 
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with initial conditions P^r = . . . = -P_i = and Pq = 1. 

There are several ways to increase the degree of multiple orthogonal poly- 
nomials, since we are working with multi-indices. Another recurrence rela- 
tion only uses the nearest neigbors of Pfi. This will be described in the next 
subsection. 

1.3 Nearest neighbor recurrence relations 

Here we briefly introduce a second recursion for multiple orthogonal poly- 
nomials. For more details we refer to pj. The nearest neighbor recurrence 
relation ([5]) connects type II multiple orthogonal polynomial Pft with the 
polynomial of degree one higher Pa+e^ and all the neighbors of degree one 
lower Ph-cj for j = 1, . . . ,r. To derive it, we proceed as follows. Since both 
Pfi and Pn+e^ ai'fi monic polynomials, the difference 

•^Pnij^^ Pfi+ekiS^') 

is a polynomial of degree < |n|. Choosing hfi^k appropriately we can also 
cancel the term containing x'"'. Hence 

xPn{x) — Pfi^gi^{x) — bfi^kPn{x) 

is a polynomial of degree < |n| — 1. It is easy to see that this polynomial 
is orthogonal to all polynomials of degree nj — 2 with respect to fij, for j = 
1, . . . r. Finally, the polynomials Pa-en ■ ■ ■ i Pn-Sr form a basis for the linear 
space of all polynomials of degree < |n| — 1 that satisfy the orthogonality 
conditions 

x^Pfi{x) dfij{x) = 0, for k < Uj — 2, and j = 1, . . . r. 
Hence we can write the polynomial xPfi{x) — Pfi+g^{x) — hnkPni 

linear 

combination of the polynomials Pa-Sii ■ ■ ■ ■, Pn-Sri with some coefficients a^j. 
This forms the so-called nearest neighbor recurrence relation 

r 

xPfi{x) = + h.n,kPn{x) + ^ anjPn-ey (5) 

There are r such relations for each = 1, . . . , r. Note that the coefficients a^^- 
do not depend on k. Indeed, the coefficients be computed directly 
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from this recursion using the orthogonahty conditions (j2]), giving 



(6) 



while the coefficients 6^^^ can be computed by muhiplying both sides of the 
equation by the multiple orthogonal polynomials of type I and applying the 
bi-orthogonality (see |1] Chapter 23, Theorem 23.1.6), giving 



Similar recursion relations hold for type I multiple orthogonal polynomi- 
als, see in [H] . 

2 Interlacing properties 

D. Kershaw proved in [5] the interlacing property for the zeros of polynomials 
orthogonal with respect to a Markov system. If the sequence 02 (a^), • • ■ 

forms an integrable Markov system on (a, b) (in particular a Chebyshev sys- 
tem) then there exists a unique polynomial g„ of degree equal to n with real, 
simple zeros in [a,b], such that 



and if Qn+i is a polynomial of degree exactly n + 1 with real, simple zeros 
satisfying the same orthogonality condition, then the zeros of qn and qn+i 
interlace. 

Later, in [2], this was generalized in the sense that the Lebesgue measure 
was replaced by an arbitrary Borel measure, which was needed to prove inter- 
lacing in the case of multiple orthogonal polynomials for a Nikishin system. 
The authors showed that if A^(cri, . . . ar) is an arbitrary Nikishin system of r 
measures, then the zeros of the type II multiple orthogonal polynomials P^j 
and Pn+Sk interlace and the zeros of and Q^j interlace, where 






g„(a;)0j(x) = 0, for z = 1, 2. . . . , n. 
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and Sj (for 1 < j < r) are measures defined by induction using 




i.e., Si = (cTi) = (Ti, S2 = {(Ti, 0-2), . . .,Sr = {di, . . . , 0"^). 

In [1] tlie authors generalized the above result, proving that the zeros of 
P^jt and PH+ee,k interlace, where the sequence of polynomials P^_fc is defined 
as the sequence of polynomials of degree rifc + . . . H-n^ with P^ o = Pn.r+i = 1, 
that is orthogonal to the measures which are given by 

dfikix) = — -— —dak{x), 

\'t^n,k-l[X)r'n,k+l[X)\ 

where Hn^k = and are defined recursively as follows, for each 



In this notation the case k = 1 corresponds to result in [2], i.e., the type II 
multiple orthogonal polynomials P^ are Ph,!- 

Using arguments as in [5l Property (P)], one can obtain the interlacing 
property for the type II multiple orthogonal polynomials with respect to mea- 
sures that form an AT system. Recall that a system of measures (/ii, . . . , fir) 
forms an AT system on [a, b] if the measures /ij are absolutely continuous 
with respect to a measure fi on [a, b], with d^j{x) = Wj{x) dfi{x), and 



is a Chebyshev system on [a,b]. The zeros of type II multiple orthogonal 
polynomials for a AT system are real and simple (e.g., [1]). The results in 
[5] and [2] can then be summarized as 

Theorem 2.1. Suppose that the measures (/ii, . . . , fir) form an AT system 
on [a, b] . Then the zeros of type II multiple orthogonal polynomials P^ and 
Pfi+ek interlace. 

Proof. Let A and B be two constants such that \A\ + |P| ^ 0. Consider 
APfi + BPfi^g^. Suppose that this polynomial has a multiple real zero at 
c e M, i.e., at least a double zero 



n we set ^nfii^x) = Ph{x) and 




for A; = 1, . . . , r. 



{Wl, XWi, . . . , X"^ "^Wl, W2 • • • , X"'' ^Wr} 



APriix) + BPrt+g^ix) = {X- C)2g(x) 
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for some polynomial Q of degree |n| — 1. Denote 



det 



Wi{xi] 
Wi(x2) 



Wi[X\ft\ 



X 



Wi[X2 



W2[Xi} 
W2(X2) 



rir — l 



Wr{Xi) 
^Wr(x2) 



for points Xi, . . . , in [a, b]. This determinant vanishes if and only if Xi = Xj 
for some We consider Wn{ function of x, and it 

changes sign only if x passes through the points Xi, . . . ,x\fi\_i. If Q is a 
polynomial of degree \n\ — 1 with real coefficients, having precisely r real 
zeros xi, . . . Xr on [a, b] at which Q changes sign then Q{x)Wr+i{x, Xi, . . . Xr) 
does not change sign on [a,b] (this is a Property (P) in [5]). In particular 



(x — c)'^Q{x)Wr+l{x, Xi, . . . , Xr) dfl{x) ^ 0. 

Using the orthogonality conditions 

x'^(x — c)'^Q{x)wj{x) dfi{x) =0, k = 0, . . .rij — 1, j = 1, 
we get 

/•b 

(x — CyQ{x)Wr+l{x, Xi, . . . , Xr) djJi{x) = 0, 



which forces a contradiction. This proves that APfi + BPfi^g^ has only simple 
zeros on the real line. Hence the linear system of equations 



Pn{x) 


Pfi+ek {x) 




'A 




0" 








B 








has only the trivial solution A = B = Q for every x G M. Therefore the 
matrix on the left hand side has non-zero determinant and by continuity and 
the behavior for large x we conclude that for all x G M 



Pn{x)Pi 



■efe 



Pl(x)P^+e-,(a:)>0. 



(9) 
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If we now fix two consecutive zeros Xk, Xk+i of P^, then P^{xk)P^{xk+i) < 0. 
From the inequality ([9]) we get 

P'fl{Xk)Pn+ek{Xk) < and P'fi{Xk+l)Pn+ek{Xk+l) < 0. 

Hence 

which means that there is at least one zero Uk of Pft+e^ between Xk and Xk+i- 
If x\ft\ is the greatest zero of Pn then P^(a;|,i|) > and from inequality Qj 
we have Prj+e'j.(x|^|) < 0. But since Pn+4 is a monic polynomial, the sign of 
Pn+e^i^) fo^^ large enough x is positive, which means that there is at least 
one zero of Pn+Sk the right side of x\fi\. Similarly there is at least one zero 
Ui of Pn+Sk on the left side of the least zero Xi of P^. Since the zeros of both 
polynomials are real and simple we have 

yj<Xj<yj+i for j = 1, . . . , |fl|. 

□ 

There is an alternative definition of an AT system: instead of zeros one 
may consider sign changes. In this case instead of using Property (P) from 
[5] in the crucial step, one uses the normality of the index n, which holds 
because (x — c)^(/ii, . . . ,/!,.) is an AT system. Then the polynomial Q of 
degree \n\ — 1 satisfies the orthogonality conditions of the multiple orthogonal 
polynomial with multi- index n for the AT-system (x — c)^(/ii, . . . , /i^), which 
is not possible because of the normality. 

In [3] the authors investigate a slightly more general notion of multiple 
orthogonal polynomials, i.e., the case of mixed type multiple orthogonal poly- 
nomials. They investigate them for two Nikishin systems. Let N{a\, . . . cr^ J 
and iV((Tf , . . . cr^j) be two Nikishin systems of ri and r2 measures respec- 
tively, which come from the same basis measure cTq = cTq. The polynomi- 
als an,o, • • • , «n,ri, where n = (rii, n2), = (ni,o, • • • , rii^n) e Z'"^+^ for 
i = 1, 2, such that 

1. the degree of Onj < riij — 1, j = 0, . . . , ri, not all identically equal to 
zero, 
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2. for i = 0, . . . , r2 they satisfy 

j i^n,o{x) + ^ a^^j{x)s\ j{x)^ dsl k{x) = 0, A; = 0, . . . , ng,/ - 1, 

where sj^^, are measures (and sj^^. its Cauchy transform) defined by 
induction using (|8]), i.e., so,o = (o"o) = o"o and Sj^^ = {(Xj, . . . , ak), 

are called mixed type multiple orthogonal polynomials. They show ([31 The- 
orem 3.5]) that the zeros of 



An,j = anfl + ^an,ii^)s]+i^i{x) and 



interlace, where = (rii + e^j, n2 + e^j), < £i < ri and < £2 < 1^2- Note 
that some special cases of this theorem are 

1. if ri = then the zeros of an,o and a^i Q interlace, i.e., in our notation, 
that the zeros of type II multiple orthogonal polynomials Pa and Pn+gk 
interlace, 

2. if r2 = and j = then the zeros of 

ri ri 

^an,i{x)sli{x) and ^ane,i{x)sli{x) 

i=0 i=0 

interlace, i.e., in our notation, that the zeros of the type I multiple 
orthogonal polynomials Qa and Qn+Sk interlace. 



2.1 Interlacing properties from the recurrence relation 

We are interested in finding the interlacing property from the nearest neigh- 
bor recurrence relations. Note that a^j 7^ in (|5]) whenever all multi-indices 
are normal and Uj > 0. Indeed these coefficients are defined as a ratio of two 
integrals Both integrals are non-zero due to normality of n and n + e^. 
If 

/ x'^'Pnix) dfij{x) = 
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then Pfi would satisfy the same orthogonahty conditions as Pn+e-i which 
shows that there exist a polynomial of degree < \n\ satisfying the orthog- 
onality conditions for a multi-index of length |ri| + 1 which contradicts the 
normality oi n + Ck- 

We assume that = if nj < for at least one j and for \n\ = we 
assume that P^ = 1. Our main result in this paper is 

Theorem 2.2. Suppose that the zeros of Pfi are real and simple, and all 
multi-indices are normal. If for all 1 < j < r and for all multi-indices n one 
has ttfi^j > whenever Uj > 0, then the zeros of Pfi and Pn+Sk interlace for 
every k with 1 < k < r. 

Proof. We will use a proof by induction on the length of n. For |n| = 1 we 
have one zero Xi of P^. Evaluating the nearest neighbor recurrence relation 
OS]) for type II multiple orthogonal polynomials at Xi, we get 



From our assumption we know that Pn-ej = 1 only for one j = jo and for 
the other j it is zero. It follows that the sign of Pn+s^ at xi is minus the sign 
of a^jjo' hence it is negative. On the other hand, since Pn+s^. is monic and of 
degree 2, the sign of Pn+e,, for large enough and small enough x is necessarily 
positive. This means that the zero of P^ is between two zeros of P^+4. 

Suppose that for |n| = m — 1 the zeros of P^ and Pfi+g^ interlace for all 
k = 1, . . . ,r. This means that Pfi+s^ at the zeros Xi of P^ has alternating 
signs, i.e.. 



Consider the case \n\ = m. Let Xi, . . . Xm be the zeros of P^. Evaluating the 
nearest neighbor recurrence relation ([5]) at Xj (1 < i < m) we get 



Denote by yi, . . . ym+i the zeros of Pft+s^- From the induction assumption we 
know that the zeros of Pfis- and P^ interlace, i.e.. 



r 





r 





10 



Therefore, since a^,- > and at least for one j with 1 < j < r we have 



Sgn Pfi+ek{Xi) 



.ly+i jf jg gven 



Both polynomials Pft and -Pj?+4 are monic and therefore for large enough 
X their sign is positive. But at the point Xm (the largest zero of P^), the 
polynomial Pn+Sk has negative sign. Therefore there is at least one zero Um+i 
of Pn+Ck such that Xm < Um+i- Similarly there is at least one zero of Pfi+Sk 
on the left hand side of the smallest zero of P,j, i.e., yi < xi. Consequently, 
since the zeros of Pn+ek are real and simple, we have exactly one zero of Pft 
between two consecutive zeros of Pn+Sk 

Vj < < Vj+i for j = 1, . . . , m. 

□ 



3 Examples 

Theorem 12.21 proves interlacing of zeros of many families of multiple orthog- 
onal polynomials, such as multiple Hermite, Charlier, Meixner of the first 
kind, Krawtchouk, and Laguerre polynomials of the second kind. However, 
there are also examples where the positivity condition (a^j- > for all n with 
rij > 0) is too strong, showing that the condition is sufficient but not neces- 
sary. For instance in the case of multiple Laguerre polynomials of the first 
kind, the coefficients not all positive. A similar situation occurs in the 

case of multiple Meixner polynomials of the second kind, since in some sense 
Meixner polynomials are the discrete analogs of the Laguerre polynomials. 
Nevertheless both examples have the interlacing property, which follows from 
the fact that the measures, with respect to which these polynomials satisfy 
orthogonality conditions, form an AT system. 

3.1 Multiple Hermite polynomials 

Multiple Hermite polynomials (see are given by the Rodrigues formula 
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where distinct real numbers. These polynomials are orthogonal 

with respect to measures (yU-i, . . . , fir) which are given by dfij{x) = ^^^^ dx 
on (— cxd,cxd). These measures form an AT system. An explicit expression 
for multiple Hermite polynomials is 



where H^^^ is the usual Hermite polynomial. From both formulas for ifS 
we get the coefficients in the nearest neighbor recurrence relation, and the 
recurrence relation is 

r 

i=i 

so that bfi^k = Ck/'^ and a^.j = nj/2. Clearly a^j- > whenever Uj > and 
therefore Theorem 12.21 gives the interlacing property for the zeros of multiple 
Hermite polynomials. 

3.2 Multiple Charlier polynomials 

Multiple Charlier polynomials (see [1]) are given by the Rodrigues formula 

r r 

where ai, . . . , > and 7^ aj whenever i ^ j, and V is the backward 
difference operator, i.e., 

V/(x) = /(x)-/(x-l). 

These polynomials are orthogonal with respect to discrete measures (/ii, . . . , /i^ 
which are given by Poisson measures 

k=0 

These measures form an AT system. An explicit expression for multiple 
Charlier polynomials is 

ni rir / \ / \ ^ 

ciw = 5: ... 5: h) . . . (;) n(-%)"'-''(-i)'''(--)w- 

ki=0 kr=0 ^ ^ j = l 
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From both formulas for C5 we get the coefficients in the nearest neighbor 
recurrence relation, and the recurrence relation is 

r 

xClix) = + (a, + \n\)Cl{x) + a,n,C|_,.(x), 

so that bfi^k = cik+\n\ and anj = CLjnj. Since aj > we have aaj > whenever 
Uj > and therefore by Theorem 12.21 we get the interlacing property for the 
zeros of multiple Charlier polynomials. 

3.3 Multiple Meixner polynomials of the first kind 

Multiple Meixner polynomials of the first kind (see are given by the 
Rodrigues formula 



r(|n| + 



r(|n| + /3)r(x + i)' 

where /3>0,0<q<1 for all i and q ^ Cj whenever i ^ j. These 
polynomials are orthogonal with respect to discrete measures (/xi, . . . , /i,.) 
which are given by 

k=0 

These measures form an AT system. We can compute an explicit expression 
for multiple Meixner polynomials of the first kind 

From both formulas for M^'^^ we can compute the coefficients in the nearest 
neighbor recurrence relation, and the recurrence relation is 

xM^/{x) = M|4(x)+(^(|n|+/3)^ + |j^jM„^"'(x) 
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so that 



These recurrence coefficients were not computed earher and appear here for 
the first time. Since < Cj < 1 and /3 > we see that a^j > whenever 
rij > and therefore by Theorem 12.21 we get the interlacing property of the 
zeros of multiple Meixner polynomials of the first kind. 



3.4 Multiple Krawtchouk polynomials 

Multiple Krawtchouk polynomials (see are type II multiple orthogonal 
polynomials that are orthogonal with respect to binomial measures 



For \n\ < N they are multiple Meixner polynomials of the first kind with 
/3 = — and q = -^^ti^ for < < 1 for all i. Hence in this case we can 
immediately write the nearest neighbor recurrence relation using [373| i.e., 



i=l 



j=i Pi ^ Vj -i 

so that 



hH,k = {N - \n\)pk + ^ - Pi), a^j = ^ ^\ 

Pj i 



i=i ^Pj-i' 



Since < < 1 and > we see that a^j- > whenever Uj > and there- 
fore Theorem 12.21 gives the interlacing property of the zeros of the Multiple 
Krawtchouk polynomials for which \n\ < N. 
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3.5 Multiple Laguerre polynomials of the second kind 

Multiple Laguerre polynomials of the second kind (see are given by the 
Rodrigues formula 

Li(x) = (-1)1-1 n n (e'^^"^^"'^^")^'"^'^^ 
i=i j=i ^ ' 

where a > — 1, Cj > and Cj 7^ Cj whenever i ^ j. These polynomials are 
orthogonal with respect to measures (/^i, . . . , Hr) which are given by dfij{x) = 
^a^-cjx |-g^ These measures form an AT system. One can compute 

an explicit expression for the multiple Laguerre polynomials of the second 
kind using Leibniz' rule: 



ki=0 kr=0 



\n\-\k\ 



1 ij=i s 



From both formulas for L'i we get the coefficients in the nearest neighbor 
recurrence relation, and the recurrence relation is 

xLUx) = L',^,Jx) + ( + ^ + + LUx) 



so that 



Jn.k 




Since a > — 1 we see that a^j- > whenever Uj > and from Theorem 12.21 
we then get the interlacing property for the zeros of multiple Laguerre of the 
second kind. 



3.6 Multiple Laguerre polynomials of the first kind 

Multiple Laguerre polynomials of the first kind (see [1]) are given by the 
Rodrigues formula 
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where ttj — Uj ^ X whenever i ^ Uj > —1. These polynomials are or- 
thogonal with respect to measures (/ii, . . . which are given by d^j{x) = 
x'^je~^dx on [0, oo). These measures form an AT system. Again one can 
compute an explicit expression for multiple Laguerre polynomials of the first 
kind using Leibniz' rule: 



X (i"i-i*i+*'+"')n;t.!(-i)i«i.«-i«i 

From both formulas for we get the coefficients in the nearest neighbor 
recurrence relation, and the recurrence relation is 

xLi[x) = L^_,,^{x) + {\n\ + l + n, + ak)mx) 

r r 



rij aj ^ 



rii + ai - rij - aj ' 
so that 

r 

ai — rij — aj 



hfi,k = |ri| + 1 + rifc + Ofc, a^j = nj{nj + Uj) J]^ 



Ui + ai- rij - aj 



Observe that j ^ and therefore we cannot apply Theorem I2.2[ but these 
polynomials have the interlacing property since the measures {/ii, . . . , fir) 
form an AT system so that Theorem 12.11 can be applied. It is worth noting 
that in [9] it was shown that J2]=i ^fij > 0. It is not clear whether such a 
condition is sufficient to prove the interlacing property. 

3.7 Multiple Mexiner polynomials of the second kind 

Multiple Mexiner polynomials of the second kind (see [1]) are given by the 
Rodrigues formula 

-rw - (^)'^'n(ft).^ii(f(^ 

j=i j=i J 



m + n,) JTix + 1] 
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where < c < 1 and /3j — f3j ^ Z whenever i j , Pj > 0. These polynomials 
are orthogonal with respect to discrete measures (yU-i, . . . , fir) which are given 
by 

(/3j)fcc\ 

k=0 

These measures form an AT system. An explicit expression for the multiple 
Meixner polynomials of the first kind is 



^ s-sC) (::) 



M "' \nrj (c- 1)1"! 

r j-1 

X Yl{f3j +X-Y^ /Ci)n,-fc,(-a;)|fc|- 
j=l i=l 

From both formulas for M^''^ we get the coefficients in the nearest neighbor 
recurrence relation, and the recurrence relation is 

xMtix) = m|;^,^(x)+(^ + K + /3,)^)m|'» 



so that 



3=^ ^ ^+3 J J 



1-c ■ " ■ 



1 — cY ~}- ni + l3i — Uj — l3j 



These recurrence coefficients were not computed earlier and appear here for 
the first time. Again we see that j ^ for all j. However the zeros of 
multiple Meixner of the second kind do interlace because of Theorem 12.11 For 
these multiple orthogonal polynomials one can also show that X]j=i ^n,j > 0. 
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